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ABSTRACT 

Isolated lepton momenta, in particular their directions are the most precisely measured 
quantities in pp collisions at LHC. This offers opportunities for multitude of precision 
measurements. It is of practical importance to verify if precision measurements with lep- 
tons in the final state require all theoretical effects evaluated simultaneously or if QED 
bremsstrahlung in the final state can be separated without unwanted precision loss. 

Results for final state bremsstrahlung in the decays of narrow resonances are obtained 
from the Feynman rules of QED in an unambiguous way and can be controlled with a very 
high precision. Also for resonances of non-negligible width, if calculations are appropriately 
performed, such separation from the remaining electroweak effects can be expected. 

Our paper is devoted to validation that final state QED bremsstrahlung can indeed be 
separated from the rest of QCD and electroweak effects, in the production and decay of 
Z and W bosons, and to estimation of the resulting systematic error. The quantitative 
discussion is based on Monte Carlo programs PHOTOS and SANC, as well as on KKMC which is 
used for benchmark results. We show, that for a large classes of W and Z boson observables 
as used at LHC, theoretical error on photonic bremsstrahlung is 0.1 or 0.2%, depending on 
the program options used. An overall theoretical error on QED final state radiation, i.e. 
taking into account missing corrections due to pair emission and interference with initial 
state radiation is estimated respectively at 0.2% or 0.3% again depending on the program 
option used. 
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1 Introduction 



Several of the most important measurements at LHC experiments, such as Higgs boson 
searches PQ, precision measurements of the W boson mass [21 [3j 0] or measurements of 
electroweak boson pair production [3] rely on a precise reconstruction of momenta for the 
final state leptons [6]. A substantial effort of the experimental community was devoted 
to optimize detector design and understand detector responses. Precision of 0.1% (even 
0.01% for lepton directions) is of no exception. For more details, see e.g. reviews [51 El E] - 

The QED effects of the final state radiation (FSR) play an important role in such ex- 
perimental studies. Final state bremsstrahlung is included in all simulation chains and 
indeed should be studied together with the detector response to leptons. It can not be sep- 
arated, because of infrared singularities of QED. Different approaches based on different 
theoretical simplifications are at present in use. At the level of the collinear approxima- 
tion, expressions for higher order FSR corrections are not only well defined but are in fact 
process independent. In general, QED calculations are process dependent but methods 
for obtaining results with 0(a 2 ) corrections and resummation of higher order effects are 
well established as well as techniques for evaluating theoretical errors. There is no need 
for introducing effective models. Situation is different if incoming or outgoing particles 
are unstable, an effort like documented in [S] is then needed. In case of Z and W de- 
cays if the narrow width approximation is used, the theoretical framework for QED FSR 
bremsstrahlung effects is unambiguous. In case of the Z decay (in fact for hard pro- 
cesses mediated by s-channel Z/j* exchange) the framework in which the QED final state 
bremsstrahlung is separated from other contributions is defined unambiguously as well. 
This was explored in e + e~ collisions at LEP and high precision solutions were proposed. 

For LEP I experiments Monte Carlo simulation programs based on exclusive expo- 
nentiation and featuring second order matrix elements were developed [9]. Considerable 
theoretical effort was invested in the reordering of the perturbative expansion. This opened 
a scheme for proper resummation of vacuum polarization diagrams into terms contributing 
to Z width and remaining vacuum polarization corrections [10J. The definition of final 
state bremsstrahlung was affected in a minimal way. The by-product of this effort was 
separation of the electroweak corrections into QED parts, corrections to the Z boson prop- 
agator (which have to be resummed to all orders) and remaining weak corrections which, 
with the proper choice of the calculation schemes, were small. 

In the case of LEP II this theoretical approach had to be reconsidered, because of 
the new W + W~ pair production processes. The gauge cancellation between diagrams of 
electroweak bosons exchanged in s and t channels, had to be carefully respected. The 
resummation of the dominant contribution of the vacuum polarization had to be restricted 
to the case of the constant width. It was not necessary however to reopen discussion on 
details of scheme for final state bremsstrahlung calculations because of the relatively small 
available statistics of VT-pair samples [TT| JT2] and thus limited interest in high precision 
calculations for QED FSR bremsstrahlung^. 

1 Consequences of resummation of parts of the electroweak effects are complex but will not be covered 
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It is important to stress that the QED final state effects in processes where leptons 
are produced through decays of W or Z/j* can be calculated and simulated with an 
essentially arbitrary high precision. They form a separate class of Feynman diagrams and 
developed already techniques should be explored at LHC, especially as the attractiveness 
of such an approach was confirmed [15] in the context of W mass measurement of CDF and 
DO collaborations. With the ever increasing precision, effects beyond photonic final state 
bremsstrahlung have to be of course considered as part of FSR effects as well, in particular 
those due to emission of extra lepton pairs. Also interference effects, such as initial-final 
state bremsstrahlung interference, has to be taken into account. 

The interference become an issue for separating QED FSR radiation from the rest of 
the electroweak corrections, at the level of cross section. That is why, the interference of 
photon emission from the initial state quarks and the final state lepton has to be discussed 
carefully. It is of practical importance to verify if the suppression of the interference due 
to boson's lifetime survives experimental cuts. If it is the case, then separation of effects 
due to QED final state bremsstrahlung and remaining parts of the electroweak and strong 
interaction corrections can be conveniently explored in the experimental studies. 

For this paper we assume, that in practical applications, all other corrections than 
QED final state interactions, that is remaining electroweak, and initial state hadronic 
interactions are expected to be measured with the properly defined observables. Thanks to 
this approach we will be able to achieve very competitive precision of theoretical predictions 
on the class of corrections directly affecting lepton momenta measurements. This represents 
a complementary approach to the one used in [16]. There, emphasis was put on the use of 
electroweak calculations together with all hadronic initial state interactions necessary for 
complete predictions for observables. 

Our paper is organized as follows. In Section 2 we describe two programs PHOTOS and 
SANC, and their theoretical base. Section 3 is devoted to tests of the first order QED calcu- 
lations. This is of importance in itself but also represent consistency checks of definitions of 
this part of electroweak corrections which will be considered at the amplitude level as QED 
final state bremsstrahlung (FSR). Definition of observables and calculation schemes used 
all over the paper are also given in this section. Section 4 is devoted to discussion of results 
for multiphoton emission and theoretical uncertainties in <fr* measurements. Section 5 is 
devoted to discussion for other than photonic bremsstrahlung effects which contribute to 
theoretical error of QED FSR simulated using PHOTOS or SANC. In particular, comments 
on emissions of pairs and interferences between photon emission from final state and other 
sources are given here. Section 6, Summary, closes the paper. 

Our discussion of theoretical error is limited to systematic error of QED FSR only, 
but is performed in the context of full event generation. Other effects, like orientation of 
spin state for the intermediate Ws or Z's, affecting input for calculation of QED FSR 

in this paper. Let us point to another theoretical constraint restricting naive resummation. It could be 
observed that even in the case of initial state QED bremsstrahlung for the process e + e~ — > ^ e ^e, previously 
performed step of resummation had to be partially revisited, because of amplitude featuring t-channcl W 
exchange |13j . For the complete second order matrix element of the two hard photon emission, the diagrams 
of charged Higgs boson exchange had to be taken into account |14) . 
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spin amplitudes are addressed in Section 5.3. Precision tag for the QED FSR calculation 
implemented in PHOTOS or SANC is finally given for a broad class of observables: first for 
the photonic bremsstrahlung and then for complete FSR corrections^]. 

2 Description of the programs 

Usually when discussing phenomenological processes at LHC in the context of the detector 
response one concentrates on description of the hard process and initial state QCD effects 
embedded e.g. in general purpose Monte Carlo generators featuring parton shower, models 
of underlying events and finally QCD NLO or NNLO corrections to the hard scattering 
process itself are taken into account, see e.g. [17] for review. 

In this paper we concentrate on the QED FSR effects. This determine choice of pro- 
grams which will be used by us for simulations. The final state interactions, consist of 
QED bremsstrahlung in decays of electroweak bosons (including cases of substantial virtu- 
alities) and to some degree on the other parts of weak corrections as well. Let us recall the 
massive effort of years 1980-2000 for establishing definition of calculation scheme at LEP 
[HJ CUB [20] where theoretically sophisticated and numerically essential separation of elec- 
troweak corrections into initial state, final state, vacuum polarization (including definition 
of the width) and interferences was established. 

In this context discussion of theoretical errors of all parts necessary for predictions is 
important as effects of new physics and their interfering genuine weak backgrounds is not 
straightforward to disentangle. This represents however further separate work on weak 
corrections. In the present paper we concentrate on final state radiation, especially on the 
final state photonic bremsstrahlung. 

For the purpose of this studies two programs featuring QED Final State Radiation for 
LHC applications will be first described and later used. We will start with presentation 
of the SANC system [21], because it features complete electroweak corrections as well. De- 
scription of PHOTOS [22, 23l[22J[25] implementing QED FSR photonic bremsstrahlung only, 
will follow. 

It might be useful to note that abbreviations LO and NLO in SANC and PHOTOS have 
somewhat different meaning. In the case of SANC, "LO" (the Leading Order) means just the 
tree-level Born cross section, while an exclusive-exponentiation-like notation is adopted in 
PHOTOS, where "LO" supposes exponentiation/resummation of the terms responsible after 
partial phase space integration, for leading logarithmic terms of photonic bremsstrahlung. 
Full coverage of multiphoton phase-space is assured. The same concerns "NLO": in SANC 
it means the one-loop approximation (Born plus 0(a) EW corrections), while in PHOTOS 
exponentiation of the 0{ot) result is assumed. 

2 In this paper we use the name photonic bremsstrahlung whenever we want to stress that only diagrams 
resulting from supplementing Born level amplitudes with photon lines are considered. The name final state 
radiation is used when we stress presence of additional pairs and final state interaction when we want to 
discuss separation with remaining parts of electroweak interactions. 
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2.1 SANC 



SANC is a computer system for Support of Analytic and Numeric calculations for experi- 
ments at Colliders [21] . It can be accessed through the Internet at http : / / sane . j inr . ru/ 
and at http : //pephsanc . cern. ch/. The SANC system is suited for calculations of one-loop 
QED, EW, and QCD radiative corrections to various SM processes. Automatized analytic 
calculations in SANC provide FORM and FORTRAN modules [26], which can be used as 
building blocks in computer codes for particular applications. 

For Drell-Yan-like processes within the SANC project there are implemented: 

— complete one-loop EW RC in CC [27J and NC [2B] processes; 

— photon induced DY processes [29]; 

- higher order photonic FSR in the collinear leading logarithmic approximation; 

- higher order photonic and pair FSR in the QED leading logarithmic approximation [30| 
EH E2]; 

— complete NLO QCD corrections [33j E3]j 

— Monte Carlo integrators [35] and event generators; 

- interface to parton showers in PYTHIA and HERWIG based on the standard Les Houches 
Accord format. 

Tuned comparisons with results of HORACE [36] and Z(W) GRADE for EW RC to CC 
and NC DY were performed within the scope of Les Houches '05, '07 and TEV4LHC '06 
workshops. A good agreement achieved in these comparisons confirms correctness of the 
implementation of the complete one-loop EW corrections in all these programs. 

An important feature of the SANC approach is the possibility to control and directly 
access different contributions to the observables being under consideration. In particular, 
SANC code allows to separate effects due to the final state radiation, the interference of 
initial and final radiation, the so called pure weak contributions, etc. 

Separation of the FSR contribution in the case of neutral current DY processes is 
straightforward, it naturally appears at the level of Feynman diagrams. But the corre- 
sponding separation in the case of the charged current DY processes is not so trivial. In 
general it is even not gauge invariant. For the sake of tuned comparison with PHOTOS, a 
special prescription for this separation^! was introduced into SANC. 

Let us consider a formal separation of the pure weak (PW) and QED contributions 
S PW and 5 QED to the total W + -»■ u + d decay width 

T PW + QED = T LO i5 PW + 5 QE D) _ (1) 

This process is described by 6 QED-like diagrams with virtual photon line and 3 other 
ones with real photon emission which together lead to the formula 

~ 2 /ll 7T 2 \ . - /II 3, Myy 



fiQED _ ® 



7T 



3 This prescription should be respected also in electroweak, non-QED FSR calculations used together 
with PHOTOS in practical applications. 
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where parameter fipw is the 't Hooft scale introduced for separation of QED and PW contri- 
butions. In order to separate the FSR QED contribution, we choose fipw = M w exp(— 
This value of the 't Hooft scale makes the total QED contribution to the W boson decay 
being equal to zero. This is in agreement with the corresponding treatment in PHOTOS, 
where by construction the effect of FSR to a process does not change the normalization of 
the cross section. 

2.2 PHOTOS 

Already in the era of data analysis of LEP experiments simulation of bremsstrahlung 
in decays of resonances and particles required specialized tools. In parallel to programs 
oriented toward highest possible overall precision for the whole processes in e + e~ collisions 
such as KKMC [9 J or KORALZ [38], programs dealing with decays only, gradually became of 
a broad use. The PHOTOS Monte Carlo was one of such applications [221 E3]. Naturally 
comparisons with these high precision generators became parts of test-beds for PHOTOS 
package. 

The principle of PHOTOS algorithm is to replace, on the basis of well defined rules, 
the decay vertex embedded in the event record such as HEPEVT [39J or HepMC [40J with 
the new one, where additional photons are added. Such solution, initially not aimed for 
high precision simulations, turned out to be very effective and precise as well. Phase 
space parameterization was carefully documented in [41J. Gradually for selected decays 
[4T| H2J H3j HI] , also exact matrix elements were implemented and could be activated in 
place of universal kernels^. Originally [22] . only single photon radiation was possible and 
approximations in the universal kernel were present even in the soft photon region. With 
time, multiphoton radiation was introduced [21] and then installation of exact first order 
matrix elements in W and Z decays became available with C++ implementation of PHOTOS 
[25] . The algorithm of PHOTOS is constructed in such a way, that the same function, but 
with different input kinematical variables, is used if the single photon emission or full 
multiphoton emission is requested. Such an arrangement enables tests in a rigorous first 
order emission environment. For multiphoton emission, the same kernel is used iteratively, 
thanks to the factorization properties. Technical checks are thus spared. Optimal solution 
for the iteration was chosen and verified with alternative calculations [15J H6] based on the 
second order matrix element. It was later extended to the multiphoton case for Z decays in 
Ref . [12] . Numerical tests of that paper, for distributions of generic kinematical observables 
pointed to the theoretical precision for the simulation of photon bremsstrahlung of the 0.1% 
level. 

When presenting numerical results from PHOTOS we always refer to its C++ version 
[25] with matrix elements for W and Z decays switched on. If the LO level is explicitly 
mentioned, matrix elements are replaced by universal kernels and algorithm as of FORTRAN 

4 Prior to introduction of the C++ interface matrix element kernels were available for our test only. They 
require more detailed information from the event record which was available from PHOTOS interface in 
FORTRAN. 



5 



version 2.14 [23], or higher is used. At present version [25] represents the up-to-date version 
of PHOTOS; it is available, with few technical updates, from LCG library|47j as well. 

3 Definitions and results of the first order calculations 

As a first step of our tests we have compared numerical results obtained from PHOTOS 
and from SANC programs in case of the single photon emission. These tests cross-check 
conventions used and numerical stability of the two calculations. They also verify the 
proper choice of parameters in PYTHIA8 generator, which is used to produce electroweak 
Born level events, on which PHOTOS is activated. We have monitored distributions for 
the following observables: pseudorapidity 7] of £~, transverse mass My of l~i>n pair and 
transverse momentum p? of £~ in the case of charged current and pseudorapidity rj of £~ , 
invariant mass M of £ + £~ pair and transverse momentum of £~ in the case of neutral 
current. 

The following set of input parameters was used: 

G„ = 1.6637 x 1(T 5 GeV' 2 , 

M w = 80.403 GeV, Y w = 2.091 GeV, 

M z = 91.1876 GeV, Y z = 2.4952 GeV, (2) 
V ud = 0.9738, V us = 0.2272, V cd = 0.2271, V cs = 0.9730, 
m e = 0.511 MeV, m M = 0.10566 GeV, 

and the following experiment motivated cuts were applied on momenta of the final state 
leptons: 

NC: \r](£ + )\ < 10, \t](£~)\ < 10, p ± (£ + ) > 15 GeV, p ± (£ + ) > 15 GeV, 

70 GeV < M(£ + £~) < 110 GeV; (3) 
CC: \r}(£-)\<10, Rl(0 > 0.1 GeV, p±(ve) > 0.1 GeV. 

We work in the running width scheme for W and Z boson propagators and fix value 
of the weak mixing angle: cos 6^ = M w /M z , sin 2 9w = 1 — cos 2 6^- The value of the 
electromagnetic coupling a is evaluated in the G^-scheme using the Fermi constant G^. 
the effective coupling is defined by ctG M = V^G^M^ sin 2 Ow/tt. 

To compute the hadronic cross section we have used CTEQ6L1 set of parton distribution 
functions with running factorization scale /i 2 = s, where s is squared total energy of the 
colliding partons in their center-of-mass system. 

Comparison is performed for muon and electron final states. For the electron final 
states and for each observable the bare and calo results are provided. In the calo case the 
four-momenta of the final electron and photon are combined into effective four-momentum 
of the electron when the separation A_R = a/ (Ar](e, r ))) 2 + (A0(e, r ))) 2 < 0.1. 

To check that the normalization of Born cross sections is properly adjusted between 
simulations using SANC and PHOTOS, we have completed tests at a sub-permille precision 
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Figure 1: Ratios for Born- level distributions in W — > eu decay. 
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Figure 2: Ratios for Born- level distributions in Z — > ee decay. 
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level. The corresponding results for electrons and ratios of differential cross sections are 
shown in Fig. [T]for CC and in Fig. [2] for NC. 

For SANC and PYTHIA+PHOTOS cases, the results for 0(a) corrections which are defined 
by 5 = (a°^ — a Born )/ a Born are presented on Figs. [3J1H1 As we can see from these figures 
agreement at the one-loop level was found to be excellent, at the level of 0.01% for both 
Z and W decays, once biases due to technical parameter separating hard and soft photon 
emission was properly tuned between the two calculations. 





Figure 3: 0(a) corrections for basic kinematical distributions from PYTHIA+PHOTOS and 
SANC in W -> ev decay. 



3.1 Dependence on technical parameters 

While performing tests we had to address well known technical problem of the so called 
"fco bias". In case of fixed order correction implemented into Monte Carlo algorithm, a 
threshold on energy for emitted photon, typically in the rest-frame of the decaying particle, 
has to be introduced. It regularizes infrared singularity. Below this threshold photons are 
simply integrated out and resulting contribution is summed up with virtual corrections to 
cancel the infrared singularity. Unfortunately k Q — > limit can not be reached, unless one 
accept working with negative event weights. The dependence on the technical parameter 
ko is however small for inclusive quantities such as the total cross section. The effect 
becomes enhanced on differential distributions near the resonance peaks. Such an effect 
can be observed in Fig. [9] for invariant mass in Z — > /i^j) decay. Generally this 
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Figure 4: 0(a) corrections for basic kinematical distributions from PYTHIA+PHOTOS and 
SANC in W ->• nv decay. 
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Figure 5: 0(a) corrections for basic kinematical distributions from PYTHIA+PHOTOS and 
SANC in W — > ev decay (calo electrons). 
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Figure 6: 0(a) corrections for basic kinematical distributions from PYTHIA+PHOTOS and 
SANC in Z — > ee decay. 
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Figure 7: C?(a) corrections for basic kinematical distributions from PYTHIA+PHOTOS and 
SANC in Z — > decay. 
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Figure 8: 0(oi) corrections for basic kinematical distributions from PYTHIA+PH0T0S and 
SANC in Z — > ee decay (calo electrons). 

dependence is nowadays of no interest as in practical application options of the programs 
with multiple photon emission should be used. This example however may be instructive 
for studies of ambiguities in implementation like in [IE] where PHOTOS is used for soft 
photon emission only, while the hard photon phase space is populated with the help of 
genuine POWHEG simulation for the final and initial state emissions simultaneously. 

4 Multiple photon emissions 

Let us now turn to the same observables, but calculated with the multiple photon emission 
option of SANC and PHOTOS, suitable for the actual comparisons with the data. One can see 
from Figs [T0in3| that agreement is a bit worse than for single photon case, but well within 
expected theoretical precision. One should stress that these results represent quantitative 
comparisons of different approximations used in the SANC and PHOTOS as well. In fact, 
the approximations used, contrary to the single photon case, are not identical. Results 
of PHOTOS represents the NLO calculations with exponentiation and resummation of the 
collinear terms of the first order photon emission matrix element. Results of SANC use the 
collinear leading logarithmic approximation which introduces by construction numerical 
dependence on the corresponding QED factorization scale fi 2 . Differences due to non- 

5 It is used however when PHOTOS is combined with POWHEG for soft photon emissions only. 



11 



1.01 
1.008 
1.006 
1.004 
1.002 — 



pp -> bjso = 7 TeV) 



SANC (e = 0.003)/PHOTOS (e = 0.003) 
SANC (e = 0.001 )/PHOTOS (e = 0.003) 
SANC (e = 0.0001 )/PHOTOS (e = 0.003) 
SANC (e = 0.00001 )/PHOTOS (e = 0.003) 




100 105 1 
M(nV) [GeV 



Figure 9: Ratio of invariant mass distributions from SANC and PHOTOS in Z decay as 
function of k = e = 2E 1>min / s/s. 



optimal choice of the scale fi 2 in SANC are below several permille for differential distributions 
and below 0.1% otherwise, see Figs. [T01IT31 Additional effort and care in estimation of the 
size of seemingly minor effects may be required, if further improvements on theoretical 
precision, beyond 0.1% are needed. 



4.1 Comparisons of PHOTOS with KKMC 

In Ref. |l2] we have demonstrated physics reasons behind very good, 0.1% level, agreement 
between PHOTOS and KKMC results for final state photonic bremsstrahlung. Still until now, 
only in case of Z/'j* intermediate state rigorous tests with second order QED matrix 
element Monte Carlo are available at present, and only in restricted condition of no hadronic 
activities in the initial state. For the reference simulations of qq — > Z/^* — > l + l^(n , -f) 
processes KKMC Monte Carlo [9] was usec|§. Also the Born level qq — > Z/j* — > l + l~ events 
from KKMC were used for simulations with PHOTOS. For both cases the monochromatic series 
of events with fixed virtuality of intermediate Z/j* state has been generated. This provides 
source of particularly valuable benchmarks as KKMC is the only program which features 
exclusive exponentiation combined with spin amplitudes for double photon emissions. As 
the numerical results of such quite extensive tests, more than 1000 figures are collected on 

6 This program, at present can not be used with simulation of the whole processes at LHC. An effort in 
this direction should be mentioned, see ref. but the corresponding results are not available for us at 
this moment. 
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Figure 10: Higher order corrections for basic kinematical distributions from PYTHIA+PHOTOS 
and SANC inW^-eu decay. 





Figure 11: Higher order corrections for basic kinematical distributions from PYTHIA+PHOTOS 
and SANC mW^fiv decay. 



13 




Figure 12: Higher order corrections for basic kinematical distributions from PYTHIA+PHOTOS 
and SANC in Z — > ee decay. 





Figure 13: Higher order corrections for basic kinematical distributions from PYTHIA+PHOTOS 
and SANC in Z — >• fM/J, decay. 
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our web page [50]: plots for Z — > ji + ji^{n^) and Z — > e + e~(n r y) are presented there. 

For all plots of collection [50] , before selection cuts are applied, events are boosted to the 
laboratory frame assuming fixed 4- vector of Z defined by its p z and r] z ; the two dimensional 
grid in (p z , rj z ) is constructed, with 7 bins in p z spanning region 0-50 GeV and 4 bins in 
rj z spanning from 0-2, for each point in the grid 40M events are simulated. Kinematical 
selection is applied on lepton and photon 4-momenta and kinematical distributions are 
constructed from accepted events. It is required that each lepton p l T > 20 GeV and < 
2.47, with the gap 1.37 < < 1.52 excluded. The gap in ^ corresponds to transition 
region between central and end-cap calorimeter in ATLAS detector, and is used here to 
somewhat arbitrary enhance possible effect of exclusive selection. Then, straightforward 
comparison between electron and muon cases is available; the only difference originating 
from leptons mass. For photon |?7 7 | < 2.37 and region 1.37 < |^ 7 | < 1.52 is excluded 
again, pJp > 15 GeV is required. In example shown in Fig. [TH angle between photon and 
a closer lepton is shown for p z = 9 GeV and rj z = 2. Results from KKMC and PHOTOS 
NL0 are compared. In both cases multiphoton emission is generated in Z/j* rest frame for 
uu — > Z/'-f* — > fi + [i~(n / -f) production process, with no initial state activity of any sort and 
virtuality of intermediate state M z +6 GeV=97.187 GeV, where M z is the Z-boson mass. 

For each choice of p z and rj z used to define Z/^j* state 4-momentum a set of three 
observables: angle between photon and closer lepton, directions of leptons and an overall 
acceptance rate is monitored in |50j. A general agreement of 0.1% can be concluded for all 
distributions. The results for LO restricted PHOTOS are collected there as well. 

For all cases one has to bear in mind that an overall normalization correction factors 
for cross section, like (1 + |~) in case of Z decay, have to be included when using PHOTOS 
package. 

4.2 The 0* observable 

Motivated by recent ATLAS publication [51] on precise observable (/>*, representing impor- 
tant improvement for the measurement of Z boson transverse momentum (p z ) at LHC, 
we have decided to devote section of this paper to discussion on the respective QED FSR 
corrections. 

The measurement of the Z boson transverse momentum (p z or 0*) offers a very sen- 
sitive way for studying dynamical effects of the strong interaction, complementary to the 
measurements of the associated production of bosons with jets. The knowledge of the p z 
distribution is crucial also to improve the modeling of the W boson production needed 
for a precise measurement of the W mass [4], in particular in the low p z region which 
dominates the cross section. The study of the low p z spectrum (p z < Mz), has also an 
important implications for the understanding of the Higgs signatures [1] as well as for the 
New Physics searches at the LHC [52] . 

The precision of the direct measurement of the spectrum at low p z at the LHC and 
at the Tevatron using the Z leptonic decay is limited by systematic uncertainties related 
to the knowledge and unfolding of the experiments resolution, in particular lepton energy 
scale [531 EI] • 
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Figure 14: The distribution of the angle in the laboratory frame between photon and the closer 
muon, as generated from KKMC and PHOTOS; selection cuts are applied, see the text. The 
intermediate state of virtuality Mz+6 GeV=97.187 GeV with the transverse momentum pp = 
9 GeV and pseudorapidity r] = 2 was created in the uu annihilation. The samples of 40M 
events were simulated, ratio of the surfaces under distributions is 0.9991. Relatively large (0.1%) 
difference to unity is typical for the larger values of pseudorapidity. For LO results (see web page 
[50j for extended results) the ratio for the surfaces, is of the same order, for electrons it is closer 
to 1. This spectrum feature plateau in rest frame of Z/'j* state, but boosted to laboratory frame 
plateau is deformed. 

In recent years, additional observables with better experimental resolution and less 
sensitive to experimental systematic uncertainties have been investigated [551 ESI |571 [58] . 
The optimal experimental observable to probe the low domain of Z/'-f* production at 
hadron colliders was found to be 4>* which is denned as 

<p; = tan(0 acop /2) ■ sin(6;), (4) 

where <fi a cop is an azimuthal opening angle between the two leptons and the angle G* is the 
scattering angle of the leptons with respect to the proton beam direction in the rest frame 
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of the dilepton system. The 0* angle is defined as 

cos(6;) = tanh ~ V+ ^j , (5) 

where rf~ and r/ + are the pseudorapiditities of the negatively and positively charged lepton, 
respectively. The variable 0* is correlated to the quantity p z /mu, where ran is the invariant 
mass of the lepton pair. It therefore probes the same physics as the transverse momentum 

and can be approximately related to it by p z ~ M z ■ <p* v • From the experimental point 
of view the variable 0* relies entirely on the angle reconstruction of the leptons in pair 
production, therefore on the tracking devices of high precision. 

We have studied the theoretical error on 0* distributions due to photonic bremsstrahlung 
effects. As in Section 14.11 comparison of results from PHOTOS and KKMC generators was 
performed and stored on web page [50J. Let us list the appropriate cuts and give example 
results. 

We have requested that for both leptons p l T > 20 GeV and < 2.4 . Distributions of 
dN( " Z d ^I +l ^ from KKMC and PHOTOS were monitored. As before, the monochromatic samples 

of qq — ¥ Z/Y — ¥ / + /~(rry) were generated for two virtualities: M z + 6 GeV (97.187 GeV) 
and M z — 4 GeV (87.187 GeV), for incoming up and down quarks. Generated events were 
boosted to the laboratory frame assuming fixed pf and r\ z of intermediate Z/'j* state. 
Again the grid of 7 bins in p z spanning region 0-50 GeV and 4 bins in r) z spanning from 
0-2 was used. 

In Fig. [T5]an example for one point of (p z , r] z ) grid is chosen and the <j)* distribution is 
shown. An agreement significantly better than 0.1% between KKMC and PHOTOS is observed. 
For other choices of flavour of incoming quarks, Z/'j* momentum and virtuality agreement 
is equally good [50] . 



4.3 Case of universal kernel 

So far, in all numerical results PHOTOS with first order matrix elements as available in 
C++ version were used both in Z and W decays. If only the universal kernel was used, as 
available in public FORTRAN PHOTOS version 2.14 or higher, the loss of precision would be 
noticeable, but the uncertainty calculated with respect of the total rate would remain at 
0.2% level, for photonic bremsstrahlung corrections to the shapes of distributions [50]. As 
in the NLO case, overall normalization factor has to be corrected for separately. 



5 Non photonic Final State Bremsstrahlung 

In this Section we concentrate on those effects which go beyond multiple photon emissions. 
They can be divided into three groups. Emission of additional pairs, the effect which 
certainly belong to final state emissions, effect of interference of initial-final state QED 
effects and finally all effects which are not directly related to final state radiation, but 
nonetheless may affect their matrix element calculations. 
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Figure 15: The distribution of the <fi^ in Z/j* — > ^^(nj) as generated from KKMC and 
PHOTOS; selection cuts are applied, see the text. The intermediate state of virtuality 97.187 GeV 
with the transverse momentum pt = 9 GeV and pseudorapidity r/ = 2 decaying into muon pair 
was created in the uu annihilation. The samples of 40M events were used, ratio of the surfaces 
under distributions is 0.9991. Relatively large (0.1%) difference to unity is typical for the larger 
value of pseudorapidity. For LO results (see web page [50J for extended results) the ratio for the 
surfaces, is of the same order. For electrons, both in LO and NLO cases, this ratio is closer to 1. 

5.1 Emission of pairs 

Emission of light fermion pairs should be included starting from the second order of QED, 
i.e. from the 0(a 2 ) corrections. There are two classes of diagrams which need to be 
taken into account. Emission of real pairs (Fig. [T6l) and the corresponding correction to 
the vertex (Fig. [TTj). These two effects cancel each other to a large degree due to the 
Kinoshita-Lee-Nauenberg theorem. The generic size of the effect can be expected to be of 
the order of higher order photonic bremsstrahlung corrections discussed so far. Moreover, 
it is well known from direct calculations in particular cases that the pair corrections are 
typically several times smaller than the photonic bremsstrahlung ones in the same order 
in a. Let us recall that careful studies of pair radiation effects were performed at LEP 
times [591 [60]. 

In the context of this paper we will estimate theoretical uncertainty due to neglection 
of pair effects in the Drell-Yan observables. 
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Figure 16: A typical example of real pair correction. 




Figure 17: A typical example of virtual pair correction. 

The SANC integrators allow to perform a quick calculation of the pair corrections within 
the leading logarithmic approximation. We apply here the formalism of electron structure 
(fragmentation) functions [30| l3Tj [32] which describe radiation in the approximation of 
collinear kinematics. 

The leading logarithmic approximation (LLA) was applied to take into account the 
corrections of the orders O (a n L n ), n = 2, 3. Let us remind that in the first order O (a) 
SANC has the complete calculation. The large logarithm L = \n(n 2 /mf) depends on the 
lepton mass mi and on the factorization scale fi. The latter is taken to be equal to the cms 
incoming parton energy (other choices are also possible). 

The pure photonic contribution to the non-singlet electron fragmentation function in 
the collinear leading logarithmic approximation reads: 

vi>(y,L) = ^-y) + ^{L-i)P {l \y) + \(^{L-i)^P^{y) 

+ l(^(L-l)) 2 p (3) (y) + 0(a^). (6) 

Analytic expressions for the relevant higher order splitting functions can be found in 
Refs. [3HI32]. 

For numerical evaluations of the splitting functions regularized by the plus prescription, 
we applied the phase space slicing as follows: 

= hm 1 5(1 - y) Uhy A + + 6(1 - y - A)^^} , (7) 



P«(y) 



i + y 2 
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where an auxiliary small parameter A is introduced. In actual computations we used 
A = 10 -4 and verified the independence of the numerical results from the variations of this 
parameter, definition of P 1 (y), i > 1 is given in Ref. [31j . 

The effect due to emission of real and virtual electron-positron pairs can be estimated 
using the non-singlet and singlet pair contributions to the LLA electron fragmentation 
functions: 



Expressions for the relevant singlet splitting functions R s and R s <g> P^ can be found in 



The differential cross section of the neutral current Drell-Yan process with FSR leading 
logarithmic corrections takes the factorized form 



d 5 *[ll(p + V -> X + l + {y lP + ) + l-{wp~)) = rf 3 a Born (p + p -+ X + l+(p+) + r (p-)) 



For consistency we expand the product of the fragmentation functions in a and take in only 
terms of the order O (a 2 L 2 ) and O (a 3 L 3 ). Terms coming from the product of the photonic 
and pair parts of the fragmentation functions are treated as a part of pair corrections. 

Numerical results for the contribution of pair corrections are presented in Figs. fT8lfT9l 
for W — > Iv decays. 

If further improvement in precision would be required, pair emission can be implemented 
e.g. into C++ PHOTOS generator [25] . 

We assume that in the experimental analysis, there is no specific implicit cut rejecting 
Z — > l + l~ff events. Our cut on p l T may reject some Z — > l + l~ff events, since with 
the real // pair emission, leptons I will have a somewhat softer energy spectrum. But 
this require relatively high energy to be carried out by the // pair, the dominant triple 
logarithmic term will thus cancel between contributions from diagrams of Figs. [16] and 
[P71 Let us stress that special care is needed in case of selection cuts sensitive to the soft, 
small virtuality and collinear to primary lepton pairs. Only this region of phase space 
may contribute significantly to pair corrections. Cuts affecting additional leptons of larger 
energies are thus of no concern. 

The direction of the leptons originating from Z decays (and therefore our 0* observable) 
may be affected by partial reconstruction of two leptons (or two hadronized quarks) of the 
extra pair. Resulting phenomena may be important only if the pair // is not collinear to 
any of the primary leptons I. Again this represents a non dominant effect, thus substantially 
below required precision goal of a permille level. 



Summarizing, details of pair corrections are still not important for precision tag at the 
level of 0.1-0.2%. 




x 



\p {2 \y) + ^P {1 \y) + \R s ®P {1 \y)-\R s {y) +0(« 4 L 4 ). (8) 



Ref. [ST]. 



xd Vl (V{y u L) + V^(y u L)) dy 2 (P"% 2 , L) + V^ h (y 2 , L)) . 



(9) 
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Figure 18: Higher order photonic and pair corrections (5 in %) for basic distributions from 
PYTHIA+PH0T0S and SANC in W \Tv decay. 





Figure 19: Higher order photonic and pair corrections (5 in %) for basic distributions from 
PYTHIA+PH0T0S and SANC in W~ e~v decay. 
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5.2 Initial-Final state bremsstrahlung interference 

The effect of QED-type interference between spin amplitudes for emission from the initial 
and final states represents important, even if numerically small, class of corrections. Even 
if separation of the initial and final state bremsstrahlung at the spin amplitude level is 
clear, large interference effects may make this separation of limited practical convenience. 
However in certain approaches interferences can be combined with QED final state effects 
in a convenient way. 

Before we will present numerical results, let us recall first some details of the discussion 
from LEP times, see Ref. [59]. The discussion was devoted to energies higher than resonance 
peak, that is why interference cancellations as explained e.g. in J6TJ [62] did not apply. In 
the case of discussed in this paper LHC observables, oriented on production of W and Z 
resonances, such suppression is nonetheless expected. The reason is of a physical origin; 
time separation between boson production and decay. However, as a consequence of the 
uncertainty principle this suppression can be broken with strong event selection cuts if 
these cuts would constrain the final state energies. 

In the studies of Drell-Yan processes at LHC one can restrict discussion of the interfer- 
ence to the first order only. On the technical level control of the QED O (a) interference 
contribution is realized in the SANC Mote Carlo integrator rather simply. The corresponding 
effect is computed by switching a respective flag in the code. 

For the W decays similar arguments related to intermediate state life-time apply. In this 
case however, size of corrections is calculation scheme dependent, i.e. depend on the way 
how diagrams of photon emission off W line are treated. Studies with SANC demonstrated 
that the interference is below 0.1% for LHC applications. They were completed not only 
for W but for Z as well, see e.g. [63] . 

Let us show in Fig. [201 as an example, corrections from interference to the (ft* observable 
discussed previously. The effect is small, below 0.1% for (f>* < 0.15 and rises to 0.5% for 
<fi* ~ 0.3. Recent measurement from ATLAS collaboration [51] of </>* observable extends 
to 1.3 but with large so far, statistical errors in range <p* > 0.2. If precision requirements 
would become more demanding, the effect should be included together with the FSR 
corrections. At present, size of the interference effect can be used to estimate the size of 
the corresponding theoretical uncertainty due to its omission. 

We can conclude that the initial-final state interference does not represent a problem for 
separating final state photonic bremsstrahlung from the remaining electroweak corrections 
in processes of W and Z production at LHC. This conclusion is justified for the precision 
of 0.1%, but it will have to be studied in more detail for more exclusive configurations, like 
e.g. larger values of <fi* distribution or for observables defined for off peak regions of lepton 
pair invariant masses. 

5.3 Relations with other electroweak and hadronic interactions 

In calculation of final state radiation matrix element, dependence on the direction of in- 
coming quarks is present. However one can see from [12 J that such effects due to e.g. initial 
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pp -> Z eV pp -> Z -> e + e 




Figure 20: IFI/FSR ratio in Z decay for 0* distribution. For <fi* > 0.2 interference effects 
become sizable. 

state interactions, affect numerical results for final state bremsstrahlung in a minor way, 
through the term which is in itself at the permille level, thus well below present precision 
requirement of 0.1%. 

For decays of narrow width states or when gauge symmetry can be used to separate 
phenomena from other parts of the interactions, there are no major difficulties to identify 
QED effects at the spin amplitude level. In the general case, QED FSR can be defined and 
its systematic error can be discussed as well. However, if e.g. contribution of diagrams 
featuring t-channel exchange of bosons complicate the separation, discussion of systematic 
errors of other parts of calculations may become scheme dependent. Our discussion on 
QED FSR corrections only will nonetheless be still useful for experimental applications. 

6 Summary 

We have addressed question of theoretical error for predictions of QED final state bremsstrahlung 
in decays of W and Z bosons, used in precision measurements at LHC experiments. 

Tests and comparisons of PHOTOS versus SANC programs for final state photonic 
bremsstrahlung were performed in realistic conditions. Hard processes and initial state 
hadronic interactions were simulated with the help of PYTHIA8 [6l] Monte Carlo program, 
for results with PHOTOS. For SANC its own set-up was used. Related differences in elec- 
troweak Born-level processes required careful tuning until agrement was established. 

We have started our discussion with technical tests and results obtained at the first 
order. Separation into QED FSR and remaining electroweak corrections have been stud- 
ied and verified at the spin amplitude level. We have checked that in PHOTOS and SANC 
numerically compatible, down to 0.01% precision level, schemes of such separation were 
defined. This agreement confirms proper installation of matrix elements and numerical 
stability of SANC and PHOTOS as well. Then the comparison was repeated after allowing 
multi-photon emission in both programs. An agreement necessary to estimate systematic 
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error in implementation of QED FSR photonic bremsstrahlung for the precision level of 
0.1% was found. This conclusion holds for the decays of intermediate states, produced 
from annihilation of light quarks, predominantly close to the W and Z resonance peaks 
but with tails of the distributions taken into account. The conclusion holds if NLO kernel 
is active in PHOTOS and for SANC multiphoton option. For PHOTOS with the L0 kernel 
theoretical precision is estimated to be 0.2%. This conclusion is limited to effects resulting 
from shapes of distributions and for the selection cuts discussed in the paper. In princi- 
ple, whenever new type of cuts is applied such comparison needs to be repeated. Effects 
on normalization have to be taken into account independently, either as part of genuine 
electroweak corrections (thanks to the proper choice of fipw in W decays), or as an simple 
overall factor, like (1 + |-) in case of the Z decay. 

We have estimated the size of the higher orders QED photonic bremsstrahlung cor- 
rections using other programs. The KKMC [9j [65] Monte Carlo program of LEP era, fea- 
turing exclusive exponentiation and second order matrix element for final state photonic 
bremsstrahlung was used for reference results. With the help of this program monochro- 
matic intermediate Z/'j* states of fixed virtuality were produced from annihilation of light 
quarks. This provided interesting test while grid of predefined values of (pr, rj) was popu- 
lated, in particular for <ft* observable. The differences versus NLO PHOTOS was found below 
0.1% (0.2% for PHOTOS kernel restricted to L0 only). 

Interferences of QED FSR with QED ISR was found to be below 0.1% for selected 
W and Z observables, as expected from the physics arguments. Separation of FSR radi- 
ation from the remaining electroweak effects is of a practical importance as it facilitate 
phenomenological work. Our calculation schemes are convenient from that point of view. 
Interference effect was found to be below required precision level. 

We estimate precision level of photonic final state corrections at 0.1%. With such 
precision tag separation of QED FSR from the rest of the process can be used for the sake 
of detector studies on final state leptons. Such detector studies represent also a well defined 
segment in comparison of theoretical predictions with the measured data. One exception 
is <p* distribution in region of large <p* > 0.15. Already at <p* ~ 0.3 interference reach 0.5%. 
In this region of phase space spin amplitudes for bremsstrahlung in initial and final state 
become gradually of comparable size. Emission of additional pairs was discussed as well 
and a size of effect was estimated at 0.1% level. 

We estimate an overall systematic error for FSR implementation in PHOTOS and SANC 
at 0.2% (0.3% for PHOTOS with L0 kernel). Further improvement of precision is possible, 
but requires more detailed discussion. Details of experimental acceptance have to be taken 
into account. 

At a margin of the discussion we entered investigation of dependence on scheme specific 
parameters such as electromagnetic factorization scale fi 2 or photon enegy threshold k$ = 
e used in fixed order simulations. This may be of some interest for further studies of 
uncertainties resulting from some choices of matching of FSR with hard process and/or 
initial state interactions and/or hard emission matrix elements. 
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